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We show by means of a perturbative weakly nonlinear analysis that the axisymmetric magnetoro- 
tational instability (MR1) of a viscous, resistive, incompressible rotating shear flow in a thin channel 
gives rise to a real Ginzburg-Landau equation for the disturbance amplitude. For small magnetic 
Prandtl number (P m ), the saturation amplitude is oc y/V-m and the resulting momentum transport 
scales as TZ~ X , where 1Z is the hydrodynamic Reynolds number. Simplifying assumptions, such as 
linear shear base flow, mathematically expedient boundary conditions and continuous spectrum of 
the vertical linear modes, are used to facilitate this analysis. The asymptotic results are shown to 
comply with numerical calculations using a spectral code. They suggest that the transport due to 
the nonlinearly developed MRI may be very small in experimental setups with P m <C 1. 

PACS numbers: Valid PACS appear here 
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The magnetorotational instability (MRI) has acquired 
growing theoretical and experimental interest in recent 
years, moving beyond just the astrophysical community. 
The linear MRI has been known for almost 50 years [HQ: 
Rayleigh stable rotating Couette flows of conducting flu- 
ids are destabilized in the presence of a vertical magnetic 
field if dil 2 jdr < (angular velocity decreasing outward). 
However, the MRI only acquired importance to astro- 
physics after the influential work of Balbus & Hawley Q , 
who demonstrated its viability to accretion disks (for any 
sufficiently weak field) with linear analysis supplemented 
by nonlinear numerical simulations. Enhanced transport 
(conceivably turbulent) is necessary for accretion to pro- 
ceed in disks found in a variety of astrophysical settings - 
from protostars to active galactic nuclei. The difficulty in 
identifying sources of such enhanced transport has been a 
limitation to proper theoretical understanding and mod- 
eling of these objects. Thus the MRI has been widely 
accepted as an attractive solution for enhanced trans- 
port (see the comprehensive reviews 0, H|) even though 
some questions on the nature of MRI-driven turbulence 
remain [6|. Because of the MRI's importance and some 
of its outstanding unresolved issues, several groups have 
recently undertaken projects to investigate the instability 
under laboratory conditions [7, 8, 9]. Additionally, nu- 
merical simulations specially designed for experimental 
setups have been conducted. 

In this Letter we perform a weakly nonlinear analysis 
of the MRI near threshold for a rotating flow in channel 
geometry, subject to idealized but mathematically expe- 
dient boundary conditions. This kind of approach is im- 
portant because the viability of this linear instability as 
the driver of turbulence and angular-momentum trans- 
port relies on understanding its nonlinear development 
and saturation. By complementing the above mentioned 



simulations, analytical methods remain useful to gain 
further physical insight. To facilitate an analytical ap- 
proach we make a number of simplifying assumptions so 
as to make the system amenable to well-known methods 
ToL 11 . liL 13 1 for the derivation of nonlinear envelope 
equations (i.e. with an amplitude weakly dependent on 
time and space). Additionally we assume a narrow chan- 
nel geometry, i.e. the small-gap limit of Taylor-Couette 
flow. 

The hydromagnetic equations in cylindrical coordi- 
nates [14j are applied to the neighborhood of a represen- 
tative point in the system, using the shearing box (SB) 
approximation (To| . The SB has been employed in nu- 
merous (analytical and numerical) studies of the MRI in 
accretion disks and are also appropriate for this model 
problem. Within the framework of the small shearing 
box (SSB)[lGi|, the base flow is steady and incompress- 
ible with constant pressure. The velocity V = U(x)y 
has a linear shear profile U(x) = —qUox, represent- 
ing an axisymmetric flow about a point rg, that rotates 
with a rate flo defined from the differential rotation law 
f2(r) oc flo(r /ro)~ q . We restrict ourselves to a con- 
stant initial magnetic field B = Bqz. Cartesian coordi- 
nates x, y, z are used to represent the radial (shear- wise) , 
azimuthal (stream-wise) and vertical directions, respec- 
tively. The base flow is axisymmetricaly disturbed by 
perturbations of the velocity, u = (u x ,u v ,u z ), magnetic 
field, b = (b x ,b y ,b z ), and total pressure, w. This re- 
sults, after non-dimensionalization, in the following set 
of non-linear equations: 
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where B = (b • V + Bq8 z ). Lengths were scaled by L 
(the channel width scale), time by Qq 1 = Qq 1 and the 
magnetic field by Bq = Bo. The quantities marked by 
tilde are thus dimensional. With this scaling, Bq = 1 
in (1) but we retain it for later convenience. The non- 
dimensional parameter C = B 2 /(Anp Q 2 L 2 ) = V\j\J 2 , 
is the Cowling number, measuring the dynamical impor- 
tance of the magnetic field. In addition the Reynolds 
number, 1Z = QqL 2 jv, and the magnetic Reynolds num- 
ber, lZ m = (IqL 2 /fj, are used, where v and fj are the 
microscopic kinematic viscosity and magnetic resistivity 
of the fluid, respectively. We introduce two additional 
non-dimensional parameters, which figure in what fol- 
lows: the Elsasser number, S = lZ m C = By / (Air p Q Q, Q fj) , 
and the magnetic Prandtl number, V m = lZ m /lZ. 

Linearization of (1), for perturbations of the form oc 
e st+ik x x-Hkzz ^ gj veg r i se ^ dispersion relation, 

aos 4 + ais 3 + a2S 2 + a^s + a± = 0, 

where all the a% = ai(k x ,k z ,V m ,S,C,q) but we shall 
write explicitly only 04: 
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k 2 T C(Ck 4 r V m + k 2 z S 2 ) 2 + K 2 S 2 Ck A T k 2 z -2qS A kt , (2) 



where the notation k 2 = 2(2 — q), k 2 = k 2 + k 2 is intro- 
duced. 

For given values of the parameters (we denote them 
collectively as II) there will be four distinct modes. The 
linear theory in various limits for this problem has been 
discussed in numerous publications (see @, @|) and we 
shall not elaborate upon them any further. Rather, we 
focus on situations where the most unstable mode (of the 
four) is marginal (critical) for some k z at a given value 
of k x — K. We fix K in order to focus on the marginal 
vertical dynamics within a channel (see below). 

Marginality to the MM (s = 0) implies the vanishing 
of the real coefficient 04 (as expressed above) and its 
derivative with respect to k z at some k z = Q: 

a 4 (k z = Q;K,IL) = 0, = Q; K, II) = 0. (3) 

ok z 

The second condition and ^ yield 
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If CK 2 > 2q the solutions are not physically meaningful, 
while the case CK 2 — 2q corresponds to the ideal MRI 
limit. 

We consider the properties of the MRI within the con- 
fines of an idealized model channel geometry with walls 
located at x = 0,w/K. We allow all quantities to be 
vertically periodic on a scale L z commensurate with in- 
teger multiples of2n/Q. As long as L z ^> l/Q, the limit 
of a vertically extended system (and thus a continuous 
spectrum of vertical modes) may be affected. 

We follow the fluid into instability by tuning the ver- 
tical background field away from the steady state, i.e. 
B — > 1 — e 2 A, where e <C 1 and A is an 0(1) con- 
trol parameter. This means that we are now in a po- 
sition to apply procedures of singular perturbation the- 
ory to this problem, by employing the multiple-scale (in 
z and t) method (e.g, [101]). It facilitates, by imposing 
suitable solvability conditions at each expansion order of 
the calculation in order to prevent a breakdown in the 
solutions, a derivation of an envelope equation for the 
unstable mode. Thus for any dependent fluid quantity 
F(x, z,t) we assume 



Q 2 S i (CQ 2 - Aq) + (K 2 + Q 2 )2Q 2 CS 2 {S 2 + n 2 ) 
{K 2 + Q 2 ) 2 CS 2 {6CV m Q 2 + k 2 ) + 2(K 2 + Q 2 ) 3 V m C 2 S 2 
5(X 2 + Q 2 ) 4 P m 2 C 3 l=0. (4) 



This equation, together with a^(k z = Q,K;H) — 0, can 
be solved for S and Q. The general expressions for 
S(K,V m> C,q) and Q(K,V m ,C,q) are lengthy but their 
asymptotic forms, for V m <C 1, to leading order, 
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y/16 Cq(2-q)K Q2 ^ R2 2q-CK 2 



F(x,z,t) = e l Fj(x, z, t). 



(6) 



The fact that the x and z components of the velocity 
and magnetic field perturbations can be derived from 
a streamfunction, ^(x, z) and magnetic flux function 
$(x, z) reduces the number of relevant dependent vari- 
ables F to four (u y and b y are the additional two) . These 
four dependent variables will also be used in the spectral 
numerical calculation (see below). 

For the lowest e order of the equations, resulting 
from substituting the expansions into the original PDE 
and collecting same order terms, we make the Ansatz 
Fi(x,z,t) — FiA(ez, e 2 t)e l ® z svuKx + c.c, where F\ is 
a constant (according to the variable in question), and 
where the envelope function A (an arbitrary constant 
amplitude in linear theory) is now allowed to have weak 
space (on scale Z = ez) and time (on scale T = e 2 t) 
dependencies. Because this system is tenth order in x- 
derivatives, a sufficient number of conditions must be 
specified at the edges. The smKx functional dependence 
of the Ansatz satisfies the ten boundary conditions 



d x u z = 0, 



(7) 



(2q-CK 2 ) 



2q + CK 2 ' 



(5) 



at x = 0, ixjK. These conditions are also subsequently 
satisfied to all orders of the expansion procedure because 
of the symmetry property of sinKx under the nonlin- 
ear operations of |T]). These specifications are a mixture 
of kinematic conditions on the edges: no-slip conditions 
on u x and u y with free-slip conditions on u z . The con- 
ditions are similarly mixed for the magnetic field: con- 
ducting conditions on b x and insulating conditions on b y . 
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Although these conditions are idealized, they are mathe- 
matically expedient, in that the coefficients of the result- 
ing envelope equation may be expressed analytically in 
the limit V m -C 1, and they serve to capture the salient 
features of the resulting dynamics (see below). 

The end result of the asymptotic procedure procedure 
is the well-known real Ginzburg-Landau Equation (GLE) 
which, for V m <C 1, is 



d T A = XA 



1 



VmC 



A\A\ 2 + Dd 2 A. 



(8) 



Here A = y/^A, A = £A, and the coefficients, for a value 
of q = 3/2 (corresponding to a Keplerian shear flow) are 



D = 6 



C 



3 5S 4 - 18S 2 - 32 + 2(S 2 + 16)VS 2 + 1 

4 S(S 2 + l)(4V52TT-3) : 
S 2 + 2 - 2y / S 2 ~TT) (S 2 + 1) 



5 3 (4V5 2r TT-3) 
35V5 2 + 1 - 65 



4<S 2 + l + y / S 2 + l 



(9a) 
(9b) 
(9c) 



where S(K) is as given in ([5]). The expressions for a 
general q are very long and involved. For iS 3> 1, (i.e. 
as one approaches the ideal MRI limit), these simplify 
to £ = 15/16, ( = 3/4, D = 3/2; and in general they 
remain O (1) quantities for all reasonable values of S. 

Before proceeding any further, it is important to note 
that equations (1) admit an integral relationship for the 
total disturbance energy. Application of the boundary 
conditions (Q yields a hydromagnetic extension of the 
Reynolds-Orr equation |l8[. Here it takes on the form: 



fffiy 

~dt~ 



= --^ J (u x u y - Cb x b y )dxdz 



TV 1 



-ClVbl 2 dxd 



(10) 



where |Vu| 2 = |Vu M | 2 for ^ = x,y,z and similarly 
for b. Ey is the total energy, per unit length in the 
azimuthal direction, of the disturbances in the domain, 
Ey = h J (u 2 + Cb 2 ) dxdz, and the first integral (multi- 
plied by dU/dx) represents the energy fed into the hy- 
drodynamical and magnetic disturbances by the back- 
ground shear. The integrand T = u x u y — Cb x b y is 
the relevant component of the total stress in the dis- 
turbance. Integrating it along the wall and averaging 
radially gives an estimate of the average angular mo- 
mentum transport rate inside the channel, i.e. J = 

{KMf%%fi' K Tdxdz. 

We can use our asymptotic solutions to express Ey and 
J in terms of the envelope function A(Z, T). To the first 
few leading orders in e (recalling that V m < 1) we find 
that 



Ey = e 2 a{S)C 2 A 2 + e 4 C 3 



A 4 



a and /3j are functions whose value is of O (1) for all 
reasonable values of S(K,q = 3/2). The two bracketed 
terms arise from the second order (in e) azimuthal ve- 
locity disturbances. Contributions to the total angular 
momentum transport (J = Jh + Jb) due to the hydrody- 
namic (Jh) and magnetic correlations (Jb) are, to leading 
order, 



Jh = ' - ; ' - 1 A 2 +0(e 3 ,e 2 V m ){lU 



Jb = 




A d + 0[e 3 ,e 2 V. 



(lib) 



The envelope A is found by solving the GLE, which 
is a well-studied system (see, e.g., [13j for a summary 
and references). It has three steady uniform solutions in 
one spatial dimension (here, the vertical): an unstable 
state A = 0, and two stable states A — ±A S , where 
A 2 = \((S)\V m C (note that ((S) is an 0(1) quantity). 
A s is also the saturation amplitude, because the system 
will develop towards it [l3j . 

Setting A — > A s in (jllal b) and in the expression for 
Ey, followed by some algebra, reveals that the angular 
momentum flux in the saturated state is, to leading or- 
der in Vra and e, J = e^CGS^GS^mCS -1 , while in the 
expression for the energy one term is independent of "P m , 

E v = e 4 C 3 C 2 (S)/? 2 (S) + 0(7> m ). 

The key results of this idealized analysis are thus, 



A s ~ y/V~C 



E^ 



En 



(12) 



for P m <C 1, where E is a constant, independent of V m . 

We have also performed numerical calculations, us- 
ing a 2-D spectral code to solve the original nonlinear 
equations (1), in the streamfunction and magnetic flux 
function formulation, near MRI threshold. The asymp- 
totic theory reflects the trends seen in these simulations. 
The code implements a Fourier-Galerkin expansion in 
64 x 64 modes in each of the four independent physi- 
cal variables, i.e. F = J2 n m F nm (i) smK n xe l ® mZ + c.c, 
where F = (-0, u y , b y ) T and where F nm (i) is the time- 
dependent amplitude of the particular Fourier-Galerkin 
mode in question (denoted by the indices n,m). 

We typically start with white noise initial conditions 
on F nm (0) at a level of 0.1 the energy of the background 
shear. Because of space limitations we display only a 
representative plot of runs made with S = 5.0, C = 0.05, 
(i.e. lZ m fixed), q = 3/2 and e 2 = 0.2 for a few succes- 
sively increasing values of 1Z (see Figure [lj . Note that 
with this value of e the fastest growing linear mode has 
a growth rate of ~ 0.065 in these units, and thus fully 
developed ideal MRI can not be expected. It is apparent 
that Ey saturates at a constant independent of 1Z (for 
large enough 1Z, i.e. small enough V m ), while J satu- 
rates at values that scale as 7?. -1 , as expected from our 
asymptotic analysis. 
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standing of this important instability, from laboratory 
experiments to numerical simulations to astrophysical 
disks. 

We thank Jeremy Goodman for his critically construc- 
tive comments on a previous version of this work and 
Michael Mond for insightful discussions. 
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FIG. 1: Ey (panel a) and J (panel b) as a function of time 
from numerical calculation. The black line is for 1Z = 5 x 10 3 , 
the red one for 1Z = 2 x 10 5 . The intermediate values (in blue) 
are successively for 1Z — 10 4 , 5 x 10 4 . The diamonds in panel 
(b) show the scaling predicted by our asymptotic analysis, 
which also predicts a constant final value of the disturbance 
energy, as apparent in panel (a). One box orbit = 2n in 
nondimensionalizcd units. 



The numerical and asymptotic solutions developed 
here show that in the saturated state the second order 
azimuthal velocity perturbation becomes dominant over 
all other quantities for V m <C 1: it has the spatial form 
six\2Kx, independent of z. It arises from the O (e 2 ) term 
in asymptotic expansion (|6|), and appears to be the pri- 
mary agent in the nonlinear saturation of the MRI in 
the channel. It acts anisotropically so as to modify the 
shear profile and results in a non-diagonal stress compo- 
nent (relevant for radial angular momentum transport), 
which behaves like ~ 1/72. 

These results should be compared to experiments and 
their numerical simulations. Thus far no experimental 
results for this geometry have been reported but corre- 
sponding numerical simulations do exist [201 ] . Our results 
on the scaling of transport (decreasing with increasing 72. 
and independent of 72 m for V m « 1) are qualitatively 
consistent with these numerical simulations but more 
careful comparisons are needed. In addition, issues of 
the effect of numerical resolution upon the resulting dy- 
namics must be considered, as is done in the turbulent 
dynamo problem 2l|, l22|, for example. Our analysis is 



complementary to that of Knobloch and Julien [ij] who 
have recently reported the results of an asymptotic MRI 
analysis for a developed state far from marginality. 

The trends predicted by this simplified model arc not 
qualitatively altered by the appearance of boundary lay- 
ers which arise when no-slip, perfectly conducting bound- 
ary conditions are enforced in the limit V m <C 1. Calcula- 
tions with these more realistic boundary conditions lead 
to amplitude equations, whose coefficients can not be, 
however, expressed analytically in any convenient way. 
A full presentation of such a calculation will appear in a 
forthcoming work. 

Further analytical investigations of the nonlinear MRI 
will contribute toward assembling a hierarchical under- 
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